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ABSTRACT 

All theories that attempt to explain the heating of the high temperature plasma 
observed in the solar corona are based on short bursts of energy. The intensities and 
velocities measured in the cores of quiescent active regions, however, can be steady 
over many hours of observation. One heating scenario that has been proposed to rec- 
oncile such observations with impulsive heating models is the "long nanoflare storm," 
where short duration heating events occur infrequently on many sub-resolutions strands; 
the emission of the strands is then averaged together to explain the observed steady 
structures. In this Letter, we examine the emission measure distribution predicted for 
such a long nanoflare storm by modeling an arcade of strands in an active region core. 
Comparisons of the computed emission measure distributions with recent observations 
indicate that that the long nanoflare storm scenario implies greater than 5 times more 
1 MK emission than is actually observed for all plausible combinations of loop lengths, 
heating rates, and abundances. We conjecture that if the plasma had "super coronal" 
abundances, the model may be able to match the observations at low temperatures. 

Subject headings: Sun: corona 



1. Introduction 



Determining the mechanism that heats the solar upper atmosphere is one of t he central goals i n 
observing the solar corona. Coronal heating theories predict short bursts of energy (|Klimchukll2006l ) ; 
these energy releases are often called "nanoflares," though the y are not nec essarily caused by small- 
scale magnetic reconnection events (i.e. Parker's nanoflares; |Parkerlll988l ). One int erpretation of 
this type of heating has b e en discussed in a series of papers by Cargill a nd Klimchuk (Car gill 1994 ; 



Cargill Klimchuk 



1997 



2004 



Klimchuk et al. 



2008 



Klimchukl |2009| ) . They describe nanoflare 



heating as a single short-lived heating event that occurs infrequently on an elemental strand in the 
corona; the density and temperature of the plasma associated with the strand then evolve without 
any additional heating events. (We adopt the definition of "strand" as a fundamental flux tube in 
the corona and "loop" as an observed coronal structure. A loop may be formed of a single strand, 
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which would imply the corona is being resolved, or many, sub-resolution strands.) The evolution 
and other properties of warm E UV loops have been shown to be consistent with the Cargill an d 
Klimchuk nanoflare model (e.g., IWarren et al.l 12002 ; IWinebarger et al.l 120031 ; IWarren et al.ll2003l ). 
In these loops, the nano flares occur within a narrow heating envelope (a so-called "short nanoflare 



storm," 
evolve. 



Klimchukl 120091 ) on a few sub-resolution strands, then the strands, and apparent loop, 



In th e core of active regions , however, the average apex intensity (IWarren et al. footpoint 
intensity (jAntiochos et al.ll2003l ). and Doppler and non-thermal velocities ([Brooks &: Warrenl l2009) 
are often found to be steady (with variations of ±20%) over several hours of observation, prompting 
some to suggest that the heating in active r egion cores is effectiv e ly steady, i.e., heating events occur 
frequently on a single s trand in the corona (|Sch 



Lundquist et al.ll2008l : I Winebarger et al. 



20a 



2011 



iiver et al.ll2004l: Warren Winebarger 



Klimchukl ([20oi) has instead suggested that 



2006 



2007 



these loops are heated by "long nanoflare storms" , i.e. the loops are bundles of many strands and 
each strand is heated at a different time. After a single heating event, the plasma along the strand 
is allowed to cool completely with no additional heating events on the strand. At any given time 
in a loop's evolution, it is composed of strands in all the different stages of the heating and cooling 
cycle so the dynamic properties of the plasma are "averaged out" in the observed loop. For this 
heating scenario, the observed loops would have some of the same properties as a frequently heated 
loop, such as steady intensity and velocity. 

The purpose of this Letter is to test whether the emission measure distribution predicted by 
the long nanoflare storm model is consistent with the observed emission measure distributions in 
active region cores. We first calculate the expected temperature and density evolution for nine 
representative strands with a range of lengths and heating rates typical to active region cores 
using a one-dimensional hydrodynamic simulation. We complete the simulations twice, once with 
a radiative loss function consistent with photospheric abundances and once with a radiative loss 
function consistent with coronal abundances. We then predict the expected emission measure 
distribution of the loop by assuming the loop is composed of a bundle of str ands that are each in a 
slightly different stage of the t emperature and plasma evolution (similar to iKlimchukl 120061 . 12009 ; 
Patsourakos Klimchukl 120091 ) . We characterize the "cool" (6.0 < Log T < Log T pea k) emission 
measure distribution of each loop by assuming it can be represented by EM ~ T a and finding a. 

For the photospheric radiative loss function, we find 1.6 < a < 2.0 for all simulated loops. 
For the coronal radiative loss function, we find 2.0 < a < 2.3 for all simulated loops. Finally, we 
sum the emission measure distributions for all loops to approximate a potential emission measure 
distribution of an active region core where many loops may be along the line-of-sight. We find 
the summed distributions have a = 1.8 for photospheric abundances and a = 2.0 for coronal 
abundances. We compare these results to two recent analyses of active region cores, which had 
3.1 < q < 3.4. We find that the long nanoflare storm scenario predicts > 5 times more 1MK 
emission than is actually observed in the core of these active regions. The long nanoflare storm 
also predicts more high temperature (> 6MK) emission than observed, but the observations are 
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not well constrained at those temperatures (jWinebarger et al.ll201ll ). Through this research, we 



have determined that the expected index, a, is sensitive to the abundances in the active region 
cores, with coronal abundances predicting larger values of a than photospheric abundances. We 
hypothesize that if the abundances in the active region core were "super-coronal" the long nanoflare 
storm simulation could reproduce the observed emission measure distribution. 



2. Long Nanoflare Storm Model 

In this section, we discuss the one-dimensional hydrodynamic simulations of nine representative 
strands. We consider three half-lengths (excluding the model chromosphere) of 25 Mm, 50 Mm and 
100 Mm. We assume the strands are semi-circular, perpendicular to the solar surface, and symmetric 
(only the half loop solution is solved). We parameterize the spatial and temporal dependence of 
the energy deposition as 

E H (s, t) = Eo + g{t)E F exp ( {§ \ . (1) 

where so designates the location of the peak of the heating function, cr s is the spatial width of the 
heating function, and Ep determines the maximum amplitude of the heating function. Because 
symmetry is assumed, there will be two heating locations, one along each loop leg. The function 
g(t) is chosen to be a simple triangular pulse, 

where 25 is the duration of the impulsive heating. The uniform background heating, Eq, consistent 
with ~ 0.5 MK, is always applied, so the density and temperature will eventually return to a 
steady-state equilibrium. 

In our simulations, the duration (25 = 100s), width (a s = 0.6 Mm), and heating location 
(so = 0.5L where L is the half-length) of the impulsive heating function is the same for each 
sim ulation. We vary the magnitu de of the heating, Ep, to achieve a specific equilibrium temperature 



(see lWinebarger &: Warrenll2004l ) of 4, 5, or 6MK based on the photospheric radiative loss function. 
The equilibrium temperature in an impulsive heating solution is the temperature at which the 
density, temperature, and loop length all agree with steady-state conditions. It is approximately the 
temperature at which the cooling switches from being dominated by conduction to being dominated 
by radiation. This is also approximately when the density of the solution peaks. The values of Ep 
are given in Column 3 of Table [TJ These same heating magnitudes are applied to the simulations 
using coronal abundances. 

For each of the lengths, equilibrium temperatures, and radiative loss functions (18 individual 
simulations), we solve the hydrodynamic loop equations with the Naval Research Laboratory Solar 
Flux Tube Model (NRL_SOLFTM). We adopt many of the same assumptions that were used in 
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previous simulations wit h this code, and we refer the reade r to earlier papers for additional details 
of this numerical code ([Mariskal 119871 ; iMariska et al.l Il989l ) . The radiative loss functions used in 
th e simulations are based on the atomic physics calculations with ionizatio n equilibrium described 
bv lMazzotta et al, a nd th e abundances of lGrevesse &: Sauvall (|1998l ) for photospheric abun- 



dances and 



Feldman et al.l (jl_992) for coronal abundances. 



We compute average temperature and density evolution by averaging over the top 50% of 
the coronal part of the strand where T e > 2 x 10 4 K. As an example, the evolution of the 25 Mm 
half-length strand with a heating magnitude, Ep, of 7.5 erg cm -3 s -1 is shown in Figure [TJ This 
example has an equilibrium temperature of 4MK. The left panel shows the average temperature 
as a function of time and the right panel shows the average density as a function of time for both 
photospheric and coronal abundances. As the nanoflare is turned on, the average temperature rises 
rapidly. Similarly, the density begins to increase as chromospheric material fills the strand. As 
the heating event ends, first the temperature and then the density begin to decrease. The peak 
density is higher and the cooling time longer for the photospheric solution compared to the coronal 
solution. 

Above we have described how the temperature and density evolution of representative strands 
were calculated. Here we discuss the method to combine these strands to determine properties 
of both individual loops and an active region core. First we assume that a loop is formed of 
m any strands, with a sing l e strand in each stage of t h e tem perature and density evolution, similar 
to iKlimchuk et al.l (120081 ): iPatsourakos Klimchukl (|2009j). In practical terms, this assumption 



implies a loop contains one strand for each time step of the simulation. For instance, in the 
photospheric solution shown in Figure[JJ the strand evolves for 2100 s before returning to equilibrium 
conditions. The code writes the solution every 10 s. The consequent loop, then, is formed of 210 
strands, with each strand having an average temperature and density associated with a specific 
time step of the simulation. 

The emission measure distribution of the loop is then determined by sorting the emission 
measures of the contributing strands into temperature bins. The emission measure for each strand is 
calculated from n^ds where n e is the average electron density of the strand and ds is the depth of the 
strand; this emission measure is placed in the bin containing the average temperature of the strand 
at that time step. The resulting emission measure distributions for each loop length, equilibrium 
temperature, and radiative loss assumptions are shown in Figure [2l The dotted vertical lines in 
Figure [2] mark the temperatures corresponding to the peak of the emission measure distributions; 
these peak temperatures are listed in Table [TJ For each length and equilibrium temperature, the 
photospheric solutions are shown in blue and the coronal solutions are shown in red. We scale the 
emission measure distributions so that the peak in the EM curve is 10 28 cm~ 5 . 

From Log T = 6.0 to the temperature of the peak of the emission measure distribution, the 
emission measure distribution can be approximated as EM ~ T a . We determine a by fitting the 
log of the emission measure distribution as a function of log of temperature with a linear regression 
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routine. These fits are shown in Figure [2] and the values of a are given in Table [TJ For photospheric 
abundances, a ranges from 1.6 to 2.0 and for coronal abundances a range from 2.0 to 2.3. In 
general, longer, hotter strands have larger as. The coronal solutions have larger values of a than 
the photospheric solutions. 

In the core of an active region, many loops are along the line-of-sight. To approximate the 
emission measure distribution for an active region core, we simply sum the emission measure dis- 
tributions for either the photospheric or coronal abundances for all loops shown in Figure [21 these 
emission measure distributions are shown in Figure We have normalized the emission measure 
distributions to 1 x 10 28 cm~ 5 for ease of comparison. We also compute the a for the summed 
solutions and we find a = 1.8 for photospheric abundances and a = 2.0 for coronal abundances. 
These a represent an emission measure weighted average a of the individual loops. This implies 
that when an active region core is formed of many loops, the emission measure distribution of the 
core can never have an a smaller or larger than the minimum and maximum a for each constituent 
loop. For comparison, we show two emission measure distributions from observations; we discuss 
these in detail in the following section. 



3. Discussion 

In this Letter, we have investigated the emission measure distribution predicted by the long 
nanoflare storm heating scenario. We calculated an array of one-dimensional hydrodynamic solu- 
tions for three different lengths, three different equilibrium temperatures, and two assumed radiative 
loss functions. Next we computed average temperature and density evolution (Figure [IJ by aver- 
aging over the upper 50% of the strand. We assumed a loop is a bundle of evolving strands with 
each strand in a different stage of the density and temperature evolution and generated an emission 
measure distribution for each loop. We have calculated the power-law indices, a, of the emission 
measure distributions for the temperature range 6.0 < Log T < Log T pea k for the simulated loops. 
These measurements are listed Table [TJ 

To approximate an active region core, where many loops are along the line-of-sight, we summed 
the emission measure distributions that are shown in Figure [2j The resulting emission measure 
distributions are presented in Figure The power law indices of the summed emission measure 
distributions are 1.8 and 2.0 for photospheric and coronal abundances, respectively. We compare 
the summed emission measure distributions, shown in Figure EJ to two recently published emission 
measure distributions o f active region cores . The power law indices o f the observed emission measure 



distributions are ~ 3.2 ([Warren et al.ll201ll ; IWinebarger et al.ll201ll ). When comparing these values 



with the observations, we find that the long nanoflare storm predicts > 5 times more 1 MK emission 
than observed. We also note that the long nanoflare storm also predicts more high temperature 
(> 6MK) emission than observed (see Figure [3]), but we do not stress this result for several 



reaso ns. First, the observations are not well constrained at high temperatures ([Winebarger et al 



2011 



Also the high-temperature evolution of the hydrodynamic simulations, which occurs early 
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in the simulation when conduction dominates the energy equation, are gre atly influenced by the 
location and duration of the heating event (see IWinebarger Warrenl (|2004l )). Finally, during this 
early stage, the plasma may be out of nonequilibrium ioniz a tion w hich would change the calculated 
emission measure distribution (see iBradshaw fe Klimchukl (|201ll )). 



One method of testing to see if the long nanoflare storm is viable is to examine line ratios 
observed in active region cores. Table [2] gives the expected EIS intensities for a few key EIS 
spectral lines for the summed coronal emission measure distribution shown in Figure [3j The last 
column of Table [2] is the ratio of the expected intensities to the expected intensity in the Ca XIV 
193.874 A line. For inst ance the expe c ted r atio between the Fe XII 195.119 A and the Ca XIV 
193.874 A lines is 13.56. IWarren et al.l ( 201ll ) determine the observed ratio is closer to 4.7. Also, 
the ratio between the Fe X 184.536 A and Ca XIV 193.874 A and Fe XV 284.160 A and Ca XIV 
193.874 A are predicted to be ~ 3.6 and ~ 86.69 respect ively; the observed ratios of these lines 
were ~ 0.9 and ~ 33.1 respectively in lWarren et al.l (|201ll ). 



We have presented emission measure distributions obtained from one-dimensional hydrody- 
namic simulations that combine plausible combinations of loops lengths, heating rates, and abun- 
dances simil ar to what is e xpected to be observed in the core of active regions. In early 0D nanoflare 
simulations, ICargilll ( 19941 ) found the slope of the emission measure curve calculated from 10 6 K to 
the peak of the emission measure to be 3.5< a < 4.8, though these results are based on simple 
assumptions of the cooling and draining of the plasma. 

The one-dimensional hydrodynamic code used in this study can not currently deal with in- 
creasing cross sectional area of loops and cannot be used to assess how the expansion affects the 
solutions. The slope of the DEMs of loops with expanding cross sections has been shown to get 
steeper with increasing expansion factor (for example, see Ivan den Oord et al.l 119971 ) and we are 
working with a new code that would be able to examine the cross sectional area in the future. 
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Fig. 1. — The evolution of an impulsively heated, 25 Mm half-length strand is shown. The average 
apex temperature (left) and density (right) are shown as a function of time for both photospheric 
(solid) and coronal (dashed) radiative loss functions. 
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Fig. 2. — The first, second and third rows show the emission measure distributions of 25, 50 and 
100 Mm half-length loops respectively. The first, second, and third columns show solutions with 
equilibrium temperatures of 4, 5, and 6MK, respectively. In each plot, the photospheric solution is 
shown as a blue line and the coronal solution is shown as a red line. The fit of the emission measure 
distributions are shown as well. The vertical dotted lines mark the temperatures corresponding to 
the peak of the emission measure distributions. 
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Fig. 3. — The average emission measure distribution from the different loops shown in Figure [2j 
For comparis o n, tw o emission measur e distr ibutions from two recent analyses are also shown 



Warren et al.l (120111 ) ; IWinebarger et al.l (120111 ) . We have normalized all the emission measure dis- 



tributions to 1 



10 28 cm" 5 . 



Table 1. Simulation Results 



Half 


Equilibrium 


Volumetric 


TpeakP 


ap 


TpeakC 


a c 


Length 


Temperature 


Heating Rate, Ef 


(MK) 




(MK) 




(Mm) 


(MK) 


(ergs cm -3 s _1 ) 










25 


4 


7.5 


3.2 


1.6 


4.0 


2.0 


25 


5 


16.1 


4.0 


1.8 


5.0 


2.0 


25 


6 


29.5 


4.0 


2.0 


6.3 


2.2 


50 


4 


7.2 


3.2 


1.7 


4.0 


2.0 


50 


5 


14.5 


4.0 


1.9 


5.0 


2.1 


50 


6 


27.1 


5.0 


1.9 


5.0 


2.3 


100 


4 


3.3 


3.2 


1.7 


4.0 


2.0 


100 


5 


5.2 


4.0 


1.9 


4.0 


2.3 


100 


6 


8.3 


4.0 


2.0 


5.0 


2.1 


SUM 






3.2 


1.7 


4.0 


2.0 



- 11 - 



Table 2. Expected EIS intensities 



EIS Spectral Line 




(ergs cm -2 s _1 sr _1 ) 


I / IcaXTV 

I V_< aiVl v 


Fe X 184.536A 


6.0 


1076.9 


3.62 


Fe XI 180.401A 


6.1 


3752.2 


12.63 


Fe XI 188.216A 


6.1 


1912.2 


6.43 


Fe XII 192.394A 


6.1 


1294.2 


4.35 


FeXII 195.119A 


6.1 


4031.3 


13.56 


S X 264.233A 


6.1 


227.2 


0.76 


Si X 258.375A 


6.1 


1133.9 


3.82 


Fe XIII 202.044A 


6.2 


1425.0 


4.80 


Fe XIII 203.826A 


6.2 


3930.6 


13.23 


Fe XIV 270.519A 


6.3 


1482.4 


4.99 


Fe XIV 264.787A 


6.3 


3246.6 


10.92 


Fe XV 284.160A 


6.3 


25762.0 


86.69 


Fe XVI 262.984A 


6.4 


1764.9 


5.94 


S XIII 256.686A 


6.4 


1684.2 


5.67 


Ca XIV 193.874A 


6.5 


297.2 


1.00 


Ca XV 200.972A 


6.6 


223.6 


0.75 


Ca XVI 208.604A 


6.7 


191.9 


0.65 


Ca XVII 192.858A 


6.7 


424.5 


1.43 



